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ABSTRACT 

A  previously  developed  iterative  procedure  is  applied  to 
obtain  numerical  solutions  of  the  von  Karman  equations  for 
rectangular  plates  subjected  to  a  uniform  normal  pressure.   On 
the  simply  supported  boundary  it  is  assumed  that  the  normal  mem- 
brane stress  and  the  tangential  membrane  displacement  vanish. 
Solutions  are  obtained  for  a  wide  range  of  values  of  the  loading 
parameter  and  the  aspect  ratio.   Boundary  layers  develop  both  as 
the  loading  parameter  and  the  aspect  ratio  increase.   The  stresses 
and  deflections  are  examined  and  compared  with  an  "asymptotic" 
solution  which  can  be  valid  only  in  the  interiors  of  long  plates. 
A  comparison  is  made  with  a  previously  obtained  approximate  solu- 
tion by  the  Ritz  method  for  the  square  plate. 


1.  Introduction 

A  thin,  elastic  rectangular  plate  is  deformed  by  a  uniform 
lateral  pressure  p  applied  normal  to  one  face.   The  stresses  and 
deflections  predicted  by  the  classical  linear  bending  theory  are 
valid  only  for  a  limited  range  of  lov;  pressures.   For  larger 
values  of  p  it  is  necessary  to  employ  nonlinear  theories  which 
allow  for  finite  deflections  and  account  for  the  interaction  of 
membrane  and  bending  effects.  Vie   employ  the  von  Karman  plate 
theory  [ij  which  is  obtained  from  a  nonlinear  elasticity  theory 
that  permits  snail  strains  and  small  but  finite  displacements. 
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We  consider  the  von  Karinan  equations  as  a  coupled  system  of 
four  second  order  partial  differential  equations.   The  four  de- 
pendent variables  are  proportional  to  the  normal  displacement  of 
the  midsurface,  the  stress  function  and  the  Laplacians  of  these 
quantities.   On  the  simply  supported  boundary  it  is  assumed  that 
the  normal  membrane  stress  and  the  tangential  membrane  displace- 
ment vanish.   The  stress  strain  and  strain  displaceraent  relations 
of  the  von  Karman  theory  then  imply  that  the  four  dependent  vari- 
ables vanish  on  the  boundary.   The  solutions  of  this  boundary 
value  problem  for  the  von  Karman  equations  depend  on  tv:o  para- 
meters, the  aspect  ratio  ^  and  a  loading  parameter  P  which  is 
proportional  to  p. 

In  this  paper,  numerical  solutions  are  obtained  by  an 
"accelerated"  iteration  method.   It  is  analogous  to  one  pre- 
viously employed  in  studies  of  the  nonlinear  axisymmetric  bending 
and  buckling  of  circular  plates  [2j  and  the  buckling  of  compressed 
rectangular  plates  [3] .   Each  iterate  is  determined  as  the  solu- 
tion of  Dirichlet's  problem  for  Poisson's  equation  on  a  rectangle. 
They  are  numerically  evaluated  by  approximating  the  Dirichlet 
problem  by  corresponding  difference  equations.   The  resulting 
system  of  linear  algebraic  equations,  which  is  of  quasi- 
tridiagonal  form,  is  solved  by  a  factoring  method  \]\\  .   The 
development  of  boundary  layers  is  studied  by  obtaining  solutions 
for  a  wide  range  of  parameters. 

The  present  results  are  compared  v;ith  approximate  solutions 
of  a  similar  boundary  value  problem  for  the  square  plate  obtained 
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by  a  Ritz  raethod  [5]  .   Oth.er  boundary  conditions  were  employed 
in  previous  investigations  of  the  nonlinear  bending  of  rectangular 
plates.   Approximate  solutions  were  obtained  by  perturbation 
methods  [6],  energy  methods  [?]  ,  infinite  series  e.g.  [6],  and 
finite  difference  methods  [9-ll] a   In  these  investigations 
relatively  low  values  of  P  are  considered  and  boundary  layer 
phenomena  are  not  studied.   The  finite  difference  procedures 
[9-ll] ,  which  in  some  aspects  are  related  to  our  method,  employ 
relatively  few  mesh  points, 

2.  Formulation 

The  rectangular  plate  of  thickness  t  occupies  the  region: 

0<X<a,  0<Y<b  and  |z|  <  t/2.   It  is  deformed  by  a  uniform 

pressure  p  applied  normal  to  one  face  and  directed  in  the  positive 

Z  direction.   The  displacements  of  the  midsurface  Z  = 0  in  the 

X,  Y  and  Z  directions  are  denoted  by  U(X,Y),  V(X,Y)  andW(X,Y), 

The  membrane  stresses  are  called  o  (X,Y) ,  a  (X,Y)  and  o   (X,Y) . 

X         y  xy 

The  following  dimensionless  variables  are  defined  by: 
x  =  Xb"-*-  ,  y  =  Yb"-*-  ,  0  <  x  <  i  =  a/b  ,  0  <  y  <  1  , 
u(x,y)  E  c2bt"2u(X,Y)  ,  v(x,y)  s  c^bt"2v(X,Y)  , 


(2.1) 


w(x,y)  E  ct'H7(X,Y)  ,  c^  =  12(1 -v^)  ,  p  =  c^(b)i|p  ^ 
Z^(x,y)  -%-(^)%(X»Y)  ,  2y(x,y)  =  %-  ( |)2oy(X  ,Y)  , 
2,y(x,y)  E^  (|)%y(X,Y)  , 
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where   v   is   Poisson's   ratio  and  E   is  Yoiing'a  modulus.      The    stress 

function,  f (x,y) ,  is  defined  in  terms  of  the  dimensionless 

membrane  stresses  Z  ,  2   and  2    by 

■A-    J  xy 

(2.2)       2^  =  f,yy  ,  2y  =  f,^^  ,  2^y  =  -f,^y  . 

A  comma  denotes  partial  differentiation  with  respect  to  the 
subscripted  variables  following  the  comma. 

The  dimensionless  membrane  stresses  are  given  in  terms  of 
the  dimensionless  displacements  u,  v  and  w  by  [l] 

12  12 

2  -  v2   =  u,  +-iw,    ,  2  -v2   =  V,  +4w,   , 
X    y     'x  2  'x  *   y    x     y  2  y  ' 


(2,3) 


2    =  r2(l+v)]~  [u,  +v,  +w,  M,    1  . 
xy    i-      -^      y   X   'x  'yJ 


The  moments  M  ,  M   and  M   and  the  corresponding  dimensionless 

moments  m  ,  m  and  m   are  given  in  terms  of  the  dimensionless 
x'   y      xy     ^ 

displacement,  w(x,y),  by 

C^b^     , 

(2.1^)     my(x,y)  =  -  £-^i.,y(x,Y)  =  w,^^+vw,^^  , 

.3u2 
m   (x,y)  5        I,  M   (X,Y)  =  w,    . 
^y        (l-v)Et^  ^y  ^y 

The  functions  .^'(x,y)  and  ^(x,y)  are  defined  by 
k2.S)  C1U,y)    =  /^   w(x,y)  ,  «()(x,y)  =  /\   f(x,y) 
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where  /\^   is   the   Laplacian  with  respect    to  x  and   y.      Then   the 
differential  equations   of   the   von  Kaman  theory  can  be  v^^ritten  ae 

/\<{,   =  K[w]     ,  Af   =   ^    > 

(2.6a) 

/\0    =   P  +  H[w,f]     ,  /\w   =    O  . 

Here   the   functionals  K  and  H  are   defined  as 


(2.6b) 


Kfw]    5  w,        -w,      w,         , 

^  J  'xy        'yy   'xx    ' 


Hfwjf]    =   f,      w,       +f,      w,        -2f,      w, 
*•        -*  yy     XX        *xx    *yy  'xy   'xj 


On  the   boundary  B  of   the   rectangular  region  D,   we   require   that 

(2.7a)  w  =  /\u  =  0    ,       'j 

.'     (x,y)    on  B  . 
(2.7b)  f  =  /\f  =  0    »         I 

The  conditions  (2.7a)  and  (2.[|)  imply  that  the  boundary  of  the 
plate  is  simply  supported.   The  boundary  conditions  (2,7b)  imply, 
using  (2.2)  and  (2.3),  that  along  each  edge  of  the  boundary  the 
normal  stress  vanishes  and  the  tangential  displacement  is  a 
constant  which  we  take  as  zero. 

The  complete  formulation  of  the  boundary  value  problem 
consists  of  the  differential  equations  (2.6)  and  the  boundary 
conditions  (2.7) .   The  solutions  depend  on  the  two  parameters 


i'  i 


, :.; .' .!•;■..■■.   .^  t/n 


■  5■i•^i  I^'b    DT,t'    [I   C'id  i/'    s  l.^;'- ; /;:' ^iU''!     ;'j'i''    \K..ir 


;  J  ■  > ::. } 


^\...*^'-  ^-X'''^::i 


■V  .1- 


vti^r:::    'j';,tii^  on:     :.v.'    ,.C   s:u'.i"e^v   'v  .;;j;5;j; :!./>:■.--;  t    -rivf   'ii?   ■:     ■  :i3!i:':T;v,ii     aJ.I"    nO 


no    \\U') 


P  and  y,     V/hen  a  solution  (w,f)  is  obtained  the  dimensionless 
stresses  and  moments  are  computed  from  (2,2)  and  (2.I4.).   The 
displacements  u  and  v  are  then  determined  from  (2.3)  using  the 
condition  that  the  tangential  displacements  vanish  on  the  boundary. 

3.  Numerical  Methods 

Approximate  solutions  of  the  boundary  value  problem  (2.6), 
(2.7)  are  obtained  by  an  accelerated  iteration  procedure  [3j . 
Thus,  starting  from  an  initial  estimate  w^°'(x,y)  of  the  solution 
for  fixed  P  and  ^,  a  sequence  of  iterates, 

[i)(^)(x,y),f(^^x,y),'.i^^^x,y),w(^-^l)(x,y)],  is  defined  by 
the  recursions, 

A^^^^  =  K[w(^)]  ,  ^(n)  =  0  on  B  , 

Af^^^  =  ^^^^  ,  f^^^  =  0  on  B  , 

(3.1)    A'^^''^  =  P  +H[w(^)  ,  f  (^)J  ,  Q(n)  =  0  on  B  , 

^-(n.l)  ^  Q(n)  ^  -(n.l)  =  q  on  B  , 

^(n-.l)  =e^(n-l)^(l_e)^(n)  ^ 

The  acceleration  parameter  ©  is  to  be  determined  so  that  the 
iterations  converge  as  rapidly  as  possible.   Each  iterate  in 
(3.1)  is  thus  determined  as  the  solution  of  a  linear  boundary 
value  problem  of  the  form, 
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(3.2)         /\bU,7)    =   c(x,y)  ,  g  =  0  on  B  , 

where  c(x,y)  is  a  known  function  on  the  rectangle  D.   The 
solution  g(x,y)  of  (3»2)  is  determined  approximately  using  a 
difference  method.   To  apply  this  method,  D  is  covered  by  a 
uniform  rectilinear  net  with  spacing  o  in  the  x  and  y  directions 
such  that  the  boundary  lines  of  D  coincide  with  net  lines.   The 
net  points  (x^,y.)  are  given  by 

(3o3a)   x^  =  i6  ,  Yj  =  j6  ,  i  =  0,1, ...,M  ,  j  =  0,1, ...,N  , 
where 
(3.3b)  6   =  A/lA   =  1/N  . 

At  each  net  point  (x.,y.)  interior  to  B,  it  is  assumed  that  the 
solution  of  (3.2),  g(x.,y.),  is  approximated  by  the  net  fxxnction 

J 

;g. ^i  which  satisfies  the  difference  equations 

(3.1|a)  A^Sij  =  C.J 

where 

(3.1+b)  A^gij  ^f^  I^ei^i,j+i-^Si-i,j-i*Si-i,j>i -^61+1,^-1 

"  ^^Si,j-^l-gi_l,j^Si,j.i-^gi,l,j)-20g.^.] 

is  the  nine  point  Laplace  difference  operator  [l2] •   The  net 
functions  ^c.  J  in  (3»i4-a)  are  obtained  by  replacing  the  second 
derivatives  occurring  in  the  right  hand  sides  of  (3.1)  by 
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centered  second  difference  approximations.   The  net  fionction 
vanishes  at  each  net  point  which  lies  on  B. 
The  coefficient  matrix  of  the  (M-l)X(N-l)  algebraic  equa- 
tions (3»l4-)  is  tridiagonal  with  respect  to  matrices  or  of  quasi- 
tridiagonal  form.   This  system  is  solved  by  factoring  the  matrix 
into  the  product  of  two  matrices,  one  lower  and  one  upper  tri- 
angular with  respect  to  matrices  and  then  successively  inverting 
the  triangular  systems  [h^  .   The  N-1,  (M-l)x(M~l)  inverses  that 
are  thus  required  are  obtained  by  Gauss  elimination  with  pivotal 
condensation.   For  a  fixed  cS,    these  inverses  are  computed  only 
once.   Hence,  only  matrix  multiplications  of  vectors  and  evalua- 
tions of  the  inhomogeneous  terms  in  (3«1)  are  required  to  deter- 
mine an  iterate, 

L).,  Computational  Methods" 

The  N-1  inverses  are  stored  in  the  internal  or  fast  access 
memory,  thus  limiting  the  size  of  the  net.'   Hence  for  a  square 
plate  we  are  essentially  limited  to  625  interior  net  points  . 
{6=   1/26)  and  for  a  rectangular  plate  with  J,  =  2  to  7i4-l  points 
(6=  1/20).  However,  a  sequence  of  test  calculations  with 
successively  finer  nets  indicates  that  they  are  adequate  for  the 
range  of  P  that  we  consider,  see  e.g.  Table  III  below,         '  ^ 

As  an  initial  iterate  we  use,  for  any  P  and  a  fixed  value 


All  computations  v;ere  performed  on  the  IBM-7090  and  7091^. 
computers  at  the  A. E.G.  Computing  and  Applied  Mathematics  Center 
of  the  Courant  Institute  of  Mathematical  Sciences. 

Finer  nets  can  be  employed  if  tape  storage  is  used.   This 
considerably  increases  the  computing  time. 
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of  Z,  the  numerical  solution  for  a  neighboring  value  of  P.  The 
first  initial  iterate  is  obtained  for  P  near  zero.  Small  pres- 
sures deform  the  plate  slightly  from  the  plane.   Therefore, 


we  take 


rij  (=0^1  =  1,2,...,  II-l,  J  =  1,2,...N-1  , 

and  the  iterations  converge  rapidly.   The  converged  solution  is 
then  used  as  an  initial  iterate  for  slightly  larger  P.   In  this 
manner  numerical  solutions  are  obtained  for  an  increasing 
sequence  of  loads. 

As  a  numerical  convergence  criterion  we  require 

(^•1)         R„  E   max,    Iw.^^^l)  -w!^.)|  <  e 
^       l<i<M-l    iJ      iJ  ' 

1<J<N-1 

where  e  >0  is  a  prescribed  "small"  number.  For  the  calculations, 
epsilons  were  in  the  range  1  xio"^  <  e  <  6.5  xlO"^  and  we  usually 
employed  the  smallest  e  consistent  with  the  accuracy  of  the 

computational  arithmetic.   It  should  be  noted  that   lim  R   =  0 

n — »oo 
is  only  a  necessary  condition  for  convergence.   When  the  itera- 
tions have  converged  for  specified  values  of  P  and  k,   the  dimen- 
sionless  stressed  and  moments  are  computed  from  difference  equiva- 
lents of  {2,2)    and  (2.1+), 

In  general,  the  number  of  iterations  necessary  to  satisfy 
(Ij-.l)  increases  as  P  increases  and  depends  on  the  value  of  © 
used  in  (3.1).   The  optimal  value  of  9,  ©    is  defined  as  that 
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value  of  ©  v/hich,  for  a  fixed  P  and  ^,  minimizes  the  number  of 

iterations  that  are  required  to  satisfy  (l+.l)  .   Estimates  of  9 

are  obtained  from  a  series  of  test  calculations  v/ith  coarse 

meshes.   The  estimated  9   is  found  to  decrease  as  P  increases. 

o 

The  number  of  iterations  needed  for  convergence  can  be 
decreased  if  the  acceleration  parameter  in  (3.1)  is  allowed  to 
vary  with  x  and  y.   To  illustrate  this,  numerical  experiments 
were  conducted  for  the  square  plate  with  P  =  100,000.   The  con- 
verged solution  for  P  =  95*000  was  employed  as  the  initial 
iterate  and  a  sequence  of  calculations  were  made  for  the  following 
choices  of  acceleration  parameters:   1)  6  =  .033  for  all  net 
points}  2)  9  =  .066  at  all  points  on  a  band  around  the  edge 
one  net  line  in  from  the  boundary,  9  =  .033  elsewhere;  3)  9  =  .066 
at  all  net  points  on  a  band  around  the  edge  consisting  of  the 
first  tv>/o  net  lines  in  from  the  boundary,  9  =  .033  elsewhere; 
i|.)  9  =  .066  at  the  four  corner  net  points  which  are  on  a  band 
one  line  in  from  the  boundary,  9  =  .033  elsewhere;  5)  6  =  .066 
at  the  same  four  corner  net  points  and  their  two  nearest  interior 
net  point  neighbors,  9  =  .033  elsewhere.   The  number  of  iterations 
needed  for  convergence  for  each  of  the  five  test  calculations  is 
summarized  in  Table  I, 

TABLE  I 


Test  Number 

1 

1 

2 

3 

h 

5 

Number  of  Iterations 

182 

161^ 

171 

156 

153 
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^.  Presentation  and  Discussion  of  Results" 

Numerical  solutions  wore  obtained  for  square  plates  and  a 
sequence  of  loads  in  the  range  0  <  P  <  100,000  and  for  rectangular 
plates  with  JL-   2   and  loads  in  the  range  0  <  p  <  6,000,   Solutions 
were  also  obtained  for  plates  with  Z  =  h,g   5>  6»  8*  10  and 
P  =  1,000  and  2,000,   The  final  mesh  widths  that  were  employed 
for  each  x  are  summarized  in  Table  II. 

TABLE  II 


i! 

Y     ■ 
J      1 

2 

h 

5 

6 

8 

10 

1/6 

1    26 

20 

16 

11+ 

12 

10 

10 

Only  one  numerical  solution  of  the  boundary  value  problem  was 
determined  for  each  P  and  -'-considered.   Each  solution  is 
symmetric  with  respect  to  the  lines  x  =  v2  and  y  =  1/2. 

To  estimate  the  range  of  validity  of  the  classical  linear 
bending  theory  '  [l3j  >  the  variation  with  P  for  "low"  loads  of  the 
normal  deflection  of  the  plate's  center  is  shown  in  Fig,  1,   The 
graphs  indicate  that  the  central  deflections  predicted  by  the 
linear  theory  [13J  deviate  from  those  of  the  nonlinear  theory, 
for  example,  for  the  square  plate,  by  more  than  10^  if  P  _>  500 
and  the  difference  increases  rapidly  i>rith  P, 


V 


•JH* 


All  calculations  employed  v  =  ,3. 


This   theory  is  obtained  from    (2.6a)    and    (2.7)    by   setting 
H[w,f]    =  K[w]    =  0   in    (2.6), 
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Representative  graphs  of  the  variation  with  x  at  y  =  1/2 
of  dimensionless  stresses  and  displacements  are  presented  in 
Fig's.  2  for  the  square  plate  and  an  increasing  sequence  of 
loads.   For  convenience',  the  results  are  shown  for  one  half  of 
the  plate.   The  normal  displacement.  Fig.  2a,  attains  its  maximum 
at  the  center  for  all  loads  considered.  The  plate  "flattens" 
slightly  at  the  center  for  the  larger  loads  and  the  flattening 
effect  increases  as  P  increases.  The  formation  of  boundary 
layers  as  P  increases  is  more  evident  in  the  graphs  shown  in 
Fig.  2b.   For  low  loads  m  (x  ,  1/2)  attains  its  maximum  at  the 
center.   For  larger  values  of  P  there  is  a  local  minimum  at  the 
center.   The  maximum,  which  then  occurs  between  the  center  and 
the  edge,  moves  towards  the  edge  as  P  increases.   A  similar 
boundary  layer  development  occurs  for  the  membrane  stress  S  , 

y 

see  Fig.  2c.   The  variations  of  the  shear  stress  2   and  the 
twisting  moment  m^  near  the  edge,  i.e.  at  y  =  1/26,  are  pre- 
sented in  Fig's.  3.   The  position  of  the  maximum  shear  stress 
moves  towards  the  corner  as  P  increases  i^rhich  indicates  the 
formation  of  boundary  layers. 

Boundary  layers  are  found  to  form  more  rapidly  and  at 
lower  loads  for  the  rectangular  plate  with  i' =  2  than  for  the 
square  plate. 

To  study  the  variations  of  the  solutions  with  /,  new 
variables  are  defined  by 


T  Syinmetry  conditions  imply  that  it  is  necessary  to  consider 
only  one  eighth  of  the  square  plate. 
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(5.1) 


X«  EX.//2  , 


w'(x'  ,  y)  =  w(x'  +  >/2,y)  . 


In  Fig's.  U,  w'(0,y)  and  w'(x'  ,  l/2)  are  shown  for  P  =  2000 
and  an  increasing  sequence  of  aspect  ratios.   The  "asyr.iptotic" 
solution,  indicated  by  the  dashed  curves  in  Pig's.  k->    is  obtained 
by  first  substituting  (5.1)  into  the  boundary  value  problem  (2.6) 
and  (2.7)  and  then  letting  / — )  co  ,   The  limit  plate  is  then 
contained  in  the  infinite  strip,  [x  |  <  co  ,   0  _<  y  <  1.   The 
asymptotic  solution,  w  (y) ,  which  is  independent  of  x',  is 
defined  as  the  solution  of  the  limiting  boundary  value  problem 
and  is  found  to  be. 


(5.2) 


w^(y)  =  |f(y^-2y2  +  i)  ,  f^(y)  =o  . 


It  does  not  satisfy  all  the  boundary  conditions  (2.7)  on 
X  =  +  -*^/2.   Hence  w.  may  be  an  approximation  to  the  solution 
of  (2.6),  (2,7)  for  "long"  plates  only  av/ay  from  the  ends.   More 
accurate  approximations  could  be  obtained  by  applying  appropriate 
boundary  layer  methods. 

A  comparison  is  given  in  Table  III  for  P  =  92,656.1  of 
the  center  deflections  w(l/2,  l/2)  obtained  from  a  four  term 
Ritz  calculation  for  the  square  plate  [5]  and  those  obtained  by 
the  present  method  for  successively  finer  meshes. 

TABLE  III 


Net   Size 

6  =   1/26 

C   =   1/12 

S  =   1/6 

U-  Term 
Ritz[5' 

w (1/2, 1/2) 

22.6699 

21.3592 

18.91+30 

18.7195 
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Similar  comparisons  can  be  made  for  other  loads.   The  maximiun 
deviation  between  the  Ritz  results  and  the  present  ones  for  the 
finer  mesh  ^  =■   1/26  is  not  at  the  center. 
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Captions  for  Fipjurea 

Fig.  1.   The  variations  with  P  of  the  nori^.al  deflections  at 
the  center  of  the  square,  i^  =  1,  and  rectangular 
plate,  £  =   2, 

Fig.  2a.   The  variation  of  w(x,l/2)  with  x  for  the  square  plate 
and  an  increasing  sequence  of  values  of  P. 

Fig.  2b.   The  variation  of  ra  (x,l/2)  with  x  for  the  square  plate 
and  an  increasing  sequence  of  values  of  P. 

Fig.  2c.   The  variation  of  Z  (x,l/2)  with  x  for  the  square  plate 
and  an  increasing  sequences  of  values  of  P. 

Fig.  3a.   The  variation  of  2   (x,l/26)  with  x  for  the  square 
plate  axid  an  increasing  sequences  of  values  of  P, 

Fig.  3b.   The  variation  of  ra   (x,l/26)  with  x  for  the  square 

xy 

plate  and  an  increasing  sequence  of  values  of  P , 

Fig.  L|.a,   The  variation  of  w'(0,y)  vixth.   y  for  an  increasing 

sequence  of  values  of  X  with  P  =  2000  and  a  conparison 
i^ith  the  asymptotic  solution  w  (jO  »  vihich  is  indicated 
by  the  dashed  curve,  see  {$,!)    and  (5.2). 

Fig.  li-b.   The  variation  of  w'(x',l/2)  with  x'  for  an  increasing 
sequence  of  values  of  i:  with  P  =  2000  and  a  conparison 
V  with  the  asymptotic  solution  which  is  indicated  by  the 
dashed  curve. 
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Figure 


w(x,l/2) 


Figure  2a 


-m„(x,i/2)«  10"^ 


Figur*  2b 


6r 


^y{x,l/2)    ,    10"^ 


Figure    2c 


-2xy(x,l/26)x  10'^ 


Figure  3a 


m^y(x,  1/26)  X  10-2 


Figure  3b 
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Chief,  Amed  Forces  Special 
lieapons  Project 

The  Pentagon,  Washington  25,  D.C. 

Attn:  Techn.  Info,  Div.  (2) 

Weapons  effects  Liv,  (1) 

Spec.  Field  Projects  (1) 

Blast  and  Shock  Br.  (1) 

Office  of  Secy,  of  the  Army 

The  Pentagon,  ashinrton  25,  D.C. 

Attn:  Amy  Library  (1) 


Chief  of  Staff,  Dept.  of  Army 

Washington  25,  D.C. 

Attn:     Dev.   Br,    {iiM)  Div,)  (l) 

Res.  3r.    (RM)  Div)  (1) 

Spec,  '..'eapons  Br,    (RStD)  (l) 

Office  of  Chief  of  Engineers 
Dept.   of  Army,  'Washington  25,   D.C. 
Attn:     ENG-IIL  Lib.  Br,  Adm, 

Ser.   Div,  (1) 

ENG-Vffi  Plan  Div,   Civ,  Wks      (1) 
EiIG-EB  Port,  Constr.    3r., 

Eng.  Div,  Mil,  Cons,(l) 
ENG-EA  Struc.  3r,  Lng. 

Div.  iiil.   Constr, 
ENG-NB  Spec.  Engr.  Br., 


ing^  R£JD  Liv. 


(1) 
(1) 


CO,  Engin.  Res.  Dev.  Lab. 

Fort  Jelvoir,  Virginia  (1) 

Office  of  Chief  of  Ordnance 

Dept  of  Army,  Washington  25,  D.C. 

Attn:  lies,  and  Mat.  Jr., 

^Ord.  RiD  Div.)  (1) 

CO,  Watertown  Arsenal 

Atertown,  Massachusetts 

Attn:   Lab.  Div,  (1) 

CO,   Frankford  Arsenal 

Bridesburg  Station 

Philadelphia  37,  Penna, 

Attn:  Lab.  Div,  (1) 

Office  of  Ordnance  Research 

2127  Myrtle  Dr.,  Duke  Station 

Durham,  North  Cai'olina 

Attn:  Div.  of  Engin,  Sci,        (1) 

CO,  Squier  Signal  Lab, 

tort  Monmo'^th,  NsJ, 

Attn:  Comp-,  and  Mat,  Br,        (l) 

Chief  of  iJaval  Operations 
Dept.  of  Navy,  l.ashington  25,  D.C, 
Attn:     Op  91  (1) 

Op  03EG  (1) 


Commandant,  Marine  Corps 
Headquart-^rs,   USMC 
U'ashington  25,  D.C. 


(1) 


Office  of  Chief  Signal  Officer 
Dept.  of  Army,  Wasliington  25,  D.C. 
Attn:   Engin.  and  Techn.  Div.     (1) 
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Chief,  Durean  of  Ships 

Dept,  of  Navy 

Washington  25,  D.C. 

Attn:   Code  106 

(1) 

Code  312 

(5) 

Code  320 

(1) 

Code  370 

(1) 

Code  375 

(1) 

Code  U20 

(1) 

Code  U21 

(1) 

Code  U23 

(2) 

Code  1x2$ 

(1) 

Code  hhO 

(1) 

Code  UU2 

(2) 

Code  Ul3 

(1) 

Code  525 

(1) 

Code  633 

(1) 

Chief,  -Buream  of  Aeronautics 

Dept.  of  Navy 

VJashington  25,  D.C, 

Attn:  AE-U 

(1) 

AV-3U 

(1) 

AD 

(1) 

AD-2 

(1) 

RS-7 

(1) 

RS-8 

(1) 

SI 

(1) 

AER-126 

(1) 

Chief,  3ureajn  of  Ordnance 

Dept.  of  Navy 

Washington  25,  D.C, 

Attn:  Ad3 

(1) 

Re 

(1) 

ReS 

(1) 

ReU 

(1) 

ReS5 

(1) 

ReSl 

(1) 

lien 

(1) 

Spec.  Proj.  Office,  3ur,  Ord, 
Dept,  of  Uavy 
Washington  25,  D.C. 
Attn:  Missile  Br. 

Chief,  Bur,  Yards  and  Docks 
Dept.  of  Navy 
VJashington  25,  D.C, 
Attn:  Code  D-202 

Code  D-202. 3 

Code  220 

Code  D-222 

Code  D-UIOC 

Code  D-liliO 

Code  D-500 


(2) 


(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 


Commanding  Officer  and  Director 

David  Taylor  Model  Jasin 

Washington  7,  D.C, 

Attn:  Code  liiO  (D 

Code  600  (1) 

Code  700  (1) 

Code  720  (D 

Code  725  (1) 

Code  731  (1) 

Code  7U0  (2) 

CO,  U.S.  Naval  Ordnance  Lab, 
l-^ite  Oak,  Maryland 
Attn:  Techn.  Library  (2) 

Techn.  Eval  Dept.         (1) 

Director,  Materials  Lab, 

N.Y.  Naval  Shipyard 

BrookljTi  1,  N.Y.  (1) 

CO,  Portsmouth  Naval  Shipyard 
Portsmouth,  New  Hampshire        (2) 

CO,  iiare  Island  Nav.  Shipyard 
Vallejo,  California  (2) 

CO  and  Director 

U.S.  Nav,  Electron.  Lab. 

San  Diego  52,  California         (1) 

Officer-in-Charge 
Nav,  Civ.  Engin,  Res, 
and  Eval.  Lab, 
U,S,  Nav.  oonstr.  Battal.  Center 
Port  Hueneme,  California 

Dir.,  Nav,  Air  Experimental  Sta, 
Nav,  Air  Mat.  Center,  Nav.  3ase 
Philadelphia  12,  Penna. 
Attn:  Materials  Lab, 
Structures  Lab, 


(2) 


(1) 
(1) 


Officer-in-Charge 

Underwater  cxplos.  Res,  Div, 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia 

Attn:  Dr.  A.H,  Keil  (2) 

CO,  U.S,  Nav,  Proving  Ground 
Dahlgren,  Virginia  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav,  Inspec.  of  Ordnance 

General  Dynamics  Corp,, 

Electr.  Boat  Div, 

Groton,  Connecticut  (1) 
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Supr,  of  Shipbuilding 

USN  and  Nav,  Inspec.  of  Ordnance 

Newport  News  Shipbuilding 

and  Dry  Dock  Co. 
Newport  News,  Virginia         (1) 


■3- 


Supr.  of  Shipbuilding 
USH  and  Nav,  Inspec,  of  Ordnance 
Ingalls  Shipbuilding  Corp< 
Pascagoula,  Mississippi 

CO,  U.S.  Nav.  Admin.  Unit 
KIT,  Cambridge  39,  liass, 

Officer-in-Charge 
Postgrad,  School  for 

i.'aval  Officers 
Webb  Inst,  of  Nav,  Arch. 
Crescent  Jeach  Rd. 
Glen  Cove,  L.I.,  N.Y. 

Supt»,  Nav.  Gun  Factory 
'w'ashington  25,  D.C. 

Comm.,  Kav.  Ordnance  Test  Sta. 
China  Lake,  California 
Attn:  Physics  Div. 

Mechanics  Div. 

CO,  Nav.  Crdnance  Test  Sta, 
Underwater  Ordnance  Div. 
3202  E.  Foothill  Jlvd. 
Pasadena  8,  California 
Attn:   Struc,  Div, 

CO  and  Director 

U.S.  Nav,  Engin,  Lxp,  Ctation 

Annapolis,  Maryland 

Supt,  U.S.  i!av.  Postgrad,  School 
Monterey,  California 

Comin.  Marine  Corps  Schools 

Quantico,  Virginia 

Attn:  Dir.,  MC  Dev,  Center 

Conm  Gen,,  USAF 
Washington  25,  D.C. 
Attn:   Res,  and  Dev.  Div, 

CO,  Air  iiaterial  Command 
VJright-Patterson  AF3,  Ohio 
Attn:  HCREX-D 

Structures  Div, 

CO,  USAF  Inst,  of  Technology 
Viright-Patterson  AF'B,  Ohio 
Attn:  Chief,  Appl,  liech. 
Group 


(1) 
(1) 


(1) 
(1) 


(1) 
(1) 


(1) 


(1) 


(1) 


(1) 


(1) 


(1) 
(1) 


(1) 


Director  of  Intelligence 
Keadquairters,  USAF 
Viashington  25,  D.C, 
Attn:   PV  Jr.  (Air  Targ,  Div) 


CO,  A.F,  Office  Sci,  Research 
Washington  25,  D.C, 
Attn:  Mechanics  Div, 


U,S,  Atomic  Energy  Commission 
l/ashinfton  25,  D.C, 
Attn:  Dir.  of  Research 

Dir.,  Nat,  Bur,  of  Standards 
Was'aington  25,  D,C, 
Attn:   Div,  of  Mechanics 
Engin,  Mech,  Sect, 
Aircraft  Structures 


(1) 


(1) 


(2) 


(1) 
(1) 
(1) 


Comm,,  U.S.  Coast  Guard 

1300  E,  St,,  iW 

i.ashington  25,  D.C, 

Attn:   Chief,  Test  and  Dev,  Div,  (1) 


U,S,  Maritime  Administration 
General  Admin,  Office  31dg, 
aUl  G  St,,  NW 
Ivashington  25,  D.C. 
Attn:  Chief,  Div,  Prelim, 
Design 


(1) 

Nat,  Aero,  and  Space  Admin. 

Langley  Research  Jenter 

Langley  Field,  Virginia 

Attn:   Structures  Div,  (2) 

Nat,  Aero,  and  Space  Admin, 

1512  H  St.,  NW 

'Washington  25,  D,C, 

Attn:   Loads  and  "true,  Div,       (2) 

Director,  Forest  Prod,  Lab, 

Madison,  Wisconsin  (1) 

Federal  Aviation  Agency 

Dept,  of  Commerce 

Washington  25,  ■~>»C, 

Attn:   Chief,  Air  Engin,  Div,     (1) 

Chief,  Air  and 

Equip,  Div,  (l) 

National  Science  Foundation 

1520  H,  St.,  mi 

iJashington,   D.C,  (1) 

General  Dynamic  Coit), 

Electr,  3oat  ^iv, 

Groton,  Connecticut  (1) 
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National  Academy  of  Sciences 

2101  Constitution  Ave,, 

Washington  25,  D.C. 

Attn:     Dir.,   Comm,   on  Ships 
Strxic.   Design 
Exec,  Secy,,    Joran,   on 
Undersea  ./arfare 

Newport  News  Shipbuilding  and 

Dry  Dock  Co. 
Newport  News,   Virginia 

Ingcills  Shipbuilding  Corp. 
Pascagoula,  ilississippi 

Prof.  Lynn   S,  3eedle 
Fritz  engineering  Lab. 
Lehigh  University 
Bethlehem,  Penna, 

Prof,  R.L.  Bisplinghoff 
Dept.  of  Aero,  Engineering 
I'lassachusetts  Inst,  of  Techn, 
Cambridge  39,  Massachusetts 

Prof.  H.  !.  Bleich 
Dept,  of  Civ.  Engineering 
Columbia  University 
New  York  27,  New  York 

Prof.  3. A,  3oley 
Dept  of  Civ.  iingineering 
Columbia  University 
New  York  27,  New  York 

Dr.  John  F.  Brahtz 
Southern  California  Labs 
Stanford  Research  Institute 
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Mission  St. 


South  Pasadena,  California 

Dr.  D.O.  Brush 
Struc.  Dept.  53-13 
Lockheed  Aircraft  Corp, 
Missile  Syst.  Div, 
Sunnyvale,  California 

Prof.  i3,  Budiansky 
Dept,  of  Hech,  ii^ngineering 
School  Appl,  Sciences 
Harvard  University 
Cambridge  36,  Massachusetts 

Prof.  Herbert  Deresiewicz 
Dept.  of  Civ,  Engineering 
Columbia  University 
632  W.  125th  St, 
New  York  27,  l',Y. 


(1) 
(1) 

(1) 
(1) 

(1) 

(1) 

(1) 

(1) 


(1) 


(1) 


(1) 


(1) 


Prof.  D.C.  Drucker,  Chmn. 

Div.  of  Engineering 

Brown  University 

Providence  12,  Rhode  Island       (1) 

Prof,  John  Duberg 

Dept.  of  Civ,  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof,  J,  Ericksen 

Mech,  .engineering  Deot. 

Johns  Hopkins  University 

Baltimore  18,  Maryland  (l) 

Prof,  A.C,  Eringen 

Dept.  of  Aero  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Prof.  \i;   Flugge 

Dept,  of  Mech.  Engineering 

Stanford,  California  (1) 

Mr.  M.  Goland,  VP  and  Dir. 

Southv;est  Research  Institute 

8500  Culebra  ud, 

San  Antonio  6,  Texas  (l) 

Prof.  J.N.  Goodier 

Dept  of  Mech.  Lngineering 

Stcinford  University 

Stanford,  California  (1) 

Prof.  L.K.  Goodman 
Engineering  experimental  Sta, 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Prof.  H.  Hetenyi 

The  Technical  Institute 

Northwestern  University 

Evanston,  Illinois  (1) 

Prof.  ?.G.  Hodge 

Dept.  of  Mechajiics 

Illinois  Inst,  of  Technology 

Chicago  16,  Illinois  (1) 

Prof.  N.J,  Hoff,  Head 

Div,  Aeronautical  Engineering 

Stanfoi'd  University 

Stanford,  California  (1) 

Prof,  V/.H,  Hoppmann,  II 

Dept,  of  Mechanics 

Rensselaer  Polytechnic  Inst. 

Troy,  New  York  (1) 
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Prof,  Bruce  G,  Johnston 

University  Of  Michigan 

Ann  Arbor,  Michigan  (1) 

Prof.  J.  iiempner 
Dept»  of  Aero,  Tngineering 
and  Appl.  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  Zt, 
Brooklyn  1,  N.Y.  (1) 

Prof,  H.L.  Langhaar 
Dept.  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Prof,  3. J,  Lazan,  Director 
Engineering  ixperlnental  ^ts. 
University  of  Minnesota 
Minneapolis  lij.,  liinnesota       (1) 

Prof,  L.H.  Lee 

Div,  of  Appl.  Mathematics 

Brown  University 

Providence  12,  Rhode  Island     (l) 

Prof,  George  H,  Lee,  Dir.  of  Ues, 

Rensselaer  Polytechnic  Inst, 

Troy,  N.Y.  (1) 

Hr,  S.  Levy 

GE  Electr,  Research  Lab, 

3198  Chestnut  St. 

Philadelphia,  Penna,  (l) 

Prof.  Paul  Lieber 

Geology  Department 

University  of  California 

Berkeley  U,  California         (l) 

Prof,  Joseph  Marin,  Head. 
Dept,  ^engineering  Mechanics 
College  of  :,ngin.  and  Arch, 
Pennsylvania  State  University 
University  Park,  Penna,        (l) 

Prof,  ic.D,  Hindlin 

Dept,  of  Civ,  Engineering 

Columbia  University 

632  w,  125th  3t, 

New  York  27,  New  York 

Prof.  Paul  M,  Naghdi 

Building  T-7 

College  of  Engineering 

University  of  California 

Serkeley  h,   California         (l) 


5  - 
Prof,  v/illiam  A,  Nash 
Dept  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida  (1) 

Prof,  M.M.  Newmark,  Head 

Dept  of  Civ,  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof,  E,  Orowan 
Dept,  of  Mech,  Engineering 
Massachusetts  Institute  of  Techn, 
Cambridge  39,  Massachusetts       (1) 

Prof,  ^vris  Phillips 

Dept,  of  Civ,  Engineering 

15  Prospect  St, 

Yale  University 

New  Haven,  Connecticut,  (1) 

Prof,  '.V,  Prager,  Ghmn, 

Phys.  ^ci.  Council 

•iro'^m   University 

Providence  12,  Rhode  Island       (1) 

Prof,  J.R.M,  Radok 

Dept,  of  Aero  engineering 

and  Appl,  Mechanics 
Polytechnic  Inst,  of  3rookljTi 
333  Jay  St. 
Brooklyn  1,  w.Y.  (1) 

Prof,  c,L,  Reiss 

Inst,  of  Mathematical  Sciences 

New  York  University 

U  VJashington  Place 

New  York  3,  N,Y,  (1) 

Prof,  h,  Reissner 
Dept,  of  Mathematics 
Massachusetts  Inst,  of  Technology 
Cambridge  39,  ..assachusetts       (1) 

Prof,  M,A.  Sadowsky 

Dept,  of  Mechanics 

Rennselaer  Polyteclinic  Inst, 

Troy,  Hew  York  (1) 

Prof,  ii.W,  Shaffer 

Dept,  of  Mech,  Engineering 

New  York  University 

University  Heights 

New  York  S2,   N.Y.  (1) 

Prof.  J,  Stallmeyer 

Dept,  of  Civ,  engineering 

University  of  Illinois 

Urban,  Illinois  (1) 
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Prof.  Eli  Sternberg 

Dept.  of  Mechanics 

Brown  University 

Providence  12,  Rhode  Island    (1) 

Prof.  T.Y,  Thomas 

Grad.  Inst.  Kath  and  iiech. 

Indiana  University 

Bloonington,  Indiana  (1) 

Prof,  G.:-".  Timoshenko 

School  of  .engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  A.S,  Velestos 

Dept  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Dr.  W,  »/enk 

Southwest  Research  Institute 

8500  Culbera  Rd, 

San  Antonio,  Texas  (1) 

Prof.  Dana  Young 

Yale  University 

New  Haven,  Connecticut        (1) 

Prof.  R.A.  Di  Taranto 

Dept.  of  Hech.  -i^ngineering 

Drexel  Institute 

32nd  and  Chestnut  Streets 

Philadelphia,  Penna.  (1) 

Mr.  H.K.  Koopman,  Secy, 
Vilelding  Res.  Council 

Engineering^:  Foundation 
29  W.  39th  St. 
New  York  18,  M.Y.  (2) 

Prof,    /alter  R.  Dsiniels 

School  of  Lngin.  and  Archit. 

Hov;ard  University 

Washington  1,  D.C.  (1) 

Comin.,  (Code  753) 

U.S.  Naval  Ordnance  Test  Sta. 

China  Lake  California 

Attn:  Techn.  Library  (1) 

Prof,  J.  1^,  Cermak 

Dept,  of  Civ.  engineering 

Colorado  State  University 

Fort  Collins,  Colorado        (1) 


6  - 

Prof.  If. J.  Hall 
Dept.  of  Civ,  Lngineering 
University  of  Illinois 
Urban,  Illinois 

Dr,  Hyman  Serbin 
Design  Integration  Dept. 
Hughes  Aircraft  Co, 
Culver  City,  California 

Commander,  V/ADD 
Wright-Patterson  AFB,  Ohio 
Attn:   WRC 

WWRl'iDS 

WWRMDD 

Commanding  Officer 

USWNOEU 

Kirtland  Air  Force  Base 

Albuquerque,  New  Mexico 

Attn:   Code  20 

(Dr.  J.N.  Brennan) 

Legislative  Reference  Service 
Library  of  Congress 
Viashington  25,  D.C, 
Attn:   Dr.  E.  '.'enk 

Dr.  A.  Itoss 

Aircraft  Nuclear  Propulsion  Dept, 
General  ..lectric  Co, 
Cincinnati  15,  Ohio 

Dr.  F,  Lane 

General  Applied  Science  Labs 
Stewart  and  herrick  Avenues 
Westbury,  L.I.,  N.Y. 

Commander  Ldward  Leonard 
Asst.  Navy  Representative 
HIT  Lincoln  Lab. 
Lexington  73,  Massachusetts 

Prof,  W,  Pohle 
Grad,  Math,  Dept. 
Adelphi  City,  L.I.,  N.Y. 

Dr,  Martin  Goldgerg 
Research  Dept, 
Grumann  Aircraft 
Bethpage,  L.I.,  .J.Y. 

Prof.  S.D,  Larmarajan 
Aero  "pace  De  ;t, 
San  Diego  State  University 
San  Diego  15,  California 
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